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Abstract 

The amplitude for production of a single photon associated with quark pair in the fusion 
of two off-shell gluons is calculated. The matrix element found is applied to the inclusive 
prompt photon hadroproduction at high energies in the framework of /cT-factorization QCD 
approach. The total and differential cross sections are calculated in both central and for- 
ward pseudo-rapidity regions. The conservative error analisys is performed. We used the 
unintegrated gluon distributions in a proton which were obtained from the full CCFM evo- 
lution equation as well as from the Kimber-Martin-Ryskin prescription. Theoretical results 
were compared with recent experimental data taken by the D0 and CDF collaborations at 
Fermilab Tevatron. Theoretical predictions for the LHC energies are given. 



PACS number(s): 12.38.-t, 13.85.-t 



1 Introduction 

The production of prompt (or direct) photons in hadron-hadron collisions at the Tevatron 
is a subject of intens studies [1-6]. Usually the photons are called "prompt" if they are 
coupled to the interacting quarks. The theoretical and experimental investigations of such 
processes have provided a direct probe of the hard subprocess dynamics, since the produced 
photons are largely insensitive to the effects of final-state hadronization. At the leading order, 
prompt photons can be produced via quark-gluon Compton scattering or quark-antiquark 
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annihilation [7] and so, the cross sections of these processes are strongly sensitive to the 
parton (quark and gluon) content of a protoE0. Very recently experimental data [6] on the 
inclusive prompt photon hadroproduction at the Tevatron have been presented by the D0 
collaboration. These data extend previous measurements to significantly higher values of 
the photon (namely to pt ~ 300 GeV at ^/s = 1960 GeV). In the present paper we 
will analyse the data [1-6] using the so-called /cT-factorization [9, 10] (or semihard [11, 12]) 
approach of QCD. 

The /cy-factorization approach is based on the familiar Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [13] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [14] gluon evolution equations. 
In this way, the large logarithmic terms proportional to In l/x are summed up to all orders 
of perturbation theory (in the leading logarithm approximation). It is in contrast with 
the popular Dokshitzer-Gribov-Lipatov-Altarelli-Parizi (DGLAP) [15] strategy where only 
large logarithmic terms proportional to In/i^ are taken into account. The basic dynamical 
quantity of the fc^-factorization approach is the so-called unintegrated (i.e. kr-dependent) 
gluon distribution A{x, k^, /i^) which determines the probability to find a gluon carrying the 
longitudinal momentum fraction x and the transverse momentum at the probing scale 
/i^. The unintegrated gluon distribution can be obtained from the analytical or numerical 
solution of BFKL or CCFM equations. Similar to DGLAP, to calculate the cross sections of 
any physical process the unintegrated gluon density A{x, kf,, /i^) has to be convoluted [9-12] 
with the relevant partonic cross section which has to be taken off mass shell (kj^-dependent). 
It is in clear contrast with the usual DGLAP scheme (so-called collinear factorization). Since 
gluons in the initial state are not on-shell and are characterized by virtual masses (propor- 
tional to their transverse momentum), it also assumes a modification of their polarization 
density matrix [11, 12]. In particular, the polarization vector of a gluon is no longer purely 
transversal, but acquires an admixture of longitudinal and time-like components. Other 
important properties of the /c^-factorization formalism are the additional contribution to 
the cross sections due to the integration over the k|. region above /i^ and the broadening of 
the transverse momentum distributions due to extra transverse momentum of the colliding 
partons. 

It is important that at present a complete theoretical description of the prompt photon 
production data at the Tevatron is a subject of special investigations (see, for example, [16] 
and references therein). Both the completeness of next-to-leading order (NLO) perturbative 
QCD calculations [17] and consistency of the available data have been the subject of intense 
discussions [18-25]. Despite that the NLO pQCD predictions agree with the recent high-p-r 
measurements [6] within uncertainties, there are still open questions. So, it was found [1-4] 
that the shape of the measured cross sections as a function of photon transverse energy E-r 
is poorly described by the NLO pQCD calculations: the observed E-r distribution is steeper 
than the predictions of perturbative QCD. These shape differences lead to a significant 
disagreement in the ratio of cross sections calculated at different center-of-mass energies 
^/s = 630 GeV and ^/s = 1800 GeV as a function of scaling variable xt = 2i?J/y^. It was 
claimed [2, 3] that the disagreement in the xt ratio is difficult to explain with conventional 
theoretical uncertainties connected with the scale dependence and parametrizations of the 
parton distributions. The origin of the disagreement has been attributed to the effect of 

^Also observed photon may arise from so called fragmentation processes [8], where a final state quark or 
gluon decays into 7. This contribution will be discussed below in Section 2. 
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initial-state soft-gluon radiation [21, 22]. In the papers [22, 26] it was shown that the 
observed discrepancy can be reduced by introducing some additional intrinsic transverse 
momentum of the incoming partons, which is usually assumed to have a Gaussian-like 
distribution. The average value of this increases from {k^) ~ 1 GeV to more than {k^) ~ 3 
GeV in hard-scattering processes as the ^/s increases from UA6 to Tevatron energies [22, 
25]. The importance of including the gluon emission through the resummation formalism 
was recognized and only recently this approach has been developed for inclusive prompt 
photon production [26-30]. 

In the framework of fc^— factorization formalism the treatment of fc^-enhancement in the 
inclusive prompt photon hadroproduction at Tevatron suggests a modification of the above 
simple fcr smearing picture. In this approach the transverse momentum of incoming par- 
tons is generated in the course of non-coUinear parton evolution under control of relevant 
evolution equations. In paper [18] the Kimber-Martin-Ryskin (KMR) formalism [31] (in the 
double leading logarithmic approximation) was applied to study the role of the both non- 
perturbative and perturbative components of partonic k^p in describing of the observed E-r 
spectrum. Note that this formalism is based on the standard DGLAP evolution. In calcu- 
lations [18] the usual on-shell matrix elements of quark-gluon fusion and quark-antiquark 
annihilation were evaluated with precise off-shell kinematics. In our previous paper [32] we 
have used the proper off-shell expressions for the matrix elements of these partonic subpro- 
cesses and also the KMR-constructed unintegrated parton densitie^ (which were evaluated 
independently from other authors). Our predictions for the inclusive prompt photon pro- 
duction agree well with available D0 [1, 2] and CDF [3-5] experimental data at Tevatron 
in both central and forward pseudo-rapidity regions. Perfect agreement was found also in 
the ratio of two cross sections calculated at -y/i = 630 GeV and -y/i = 1800 GeV. This ratio 
shows specific effect connected with off-shell gluons in the /cT^-factorization approach. 

We note, however, that an important component of all above calculations [18, 32] is 
the unintegrated quark distribution in a proton. At present these densities are available in 
the framework of KMR approach only. It is because there are some theoretical difficulties 
in obtaining the quark distributions immediately from CCFM or BFKL equation^ (see, 
for example, review [35] for more details). As a result the dependence of the calculated 
cross sections on the non-coUinear evolution scheme has not been investigated yet. This 
dependence in general can be significant and it is a special subject of study in the /cr- 
factorization formalism. Therefore in the present paper we will try a different and more 
systematic way. Instead of using the unintegrated quark distributions and the corresponding 
quark-gluon fusion and quark-antiquark annihilation cross sections we will calculate off- 
shell matrix element of the g*g* —>■ qq'j subprocess and then will operate in terms of the 
unintegrated gluon densities only. In this way the different non-coUinear evolution schemes 
can be applied. But we note that this matrix element covers only sea quark contribution 
from the quark-gluon fusion and quark-antiquark annihilation and therefore in this case 
contribution from the valence quarks is not taken into account. However, this contribution 
is significant only at large x and therefore can be safely accounted for in the collinear LO 
approximation as additional one. 

So, in the present paper we will investigate the different prompt photon production rates 

similar scenario has been used also in [33]. 
^Unintegrated quark density was considered recently in [34]. 
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using the off-mass shell matrix elements of the g*g* qq^ subprocess and make a systematic 
comparison of our predictions to the available D0 and CDF data [1-6]. Our special goal is to 
study the sensitivity of the calculated cross sections to the non-coUinear evolution scheme. In 
the numerical calculations we will test the different sets of unintegrated gluon distributions 
which are obtained from the full CCFM equation as well as from the Kimber-Martin-Ryskin 
approach. Additionally we give some predictions for LHC conditions. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the /cT-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Finally, in Section 4, we 
give some conclusions. 



2 Theoretical framework 



2.1 Kinematics 

We start from the kinematics (see Fig. 1). Let p^^^ and p^^-* be the four-momenta of 
the incoming protons and p be the four-momentum of the produced photon. The initial 
off-shell gluons have the four-momenta ki and k2 and the final quark and antiquark have the 
four- momenta pi and p2 and masses mg, respectively. In the pp center-of-mass frame we can 
write 

p« = v^/2(l,0,0,l), p(2)=v^/2 (1,0, 0,-1), (1) 

where ^/s is the total energy of the process under consideration and we neglect the masses 
of the incoming protons. The initial gluon four-momenta in high energy limit can be written 
as 

fci = Xip^^^ + /CiT, k2 = X2P^^^ + k2T, (2) 

where kiT and k2T are their transverse four-momenta. It is important that k^^^ = —kfrp ^ 
and k|j, = — /^It 0- Ftom the conservation laws we can easily obtain the following 
conditions: 

kiT + = Pit + P2T + Pt, 

Xl^/s ^ miTe^^ + m2Te^^ + IptIc^, (3) 

X2VS = miTe~^^ + m2Te~^^ + |pT|e~^, 

where y is the rapidity of produced photon, pix and p2T are the transverse four-momenta 
of final quark and antiquark, yi, y2, uiit and m2T are their center-of-mass rapidities and 
transverse masses, i.e. = -|- pfy 



2.2 Off-shell amplitude of the g*g* qq-f subprocess 

There are eight Feynman diagrams (see Fig. 2) which describe the partonic subprocess 
9*9* ~^ QQl at the leading order in a.^ and a. Let ei, €2 and e be the initial gluon and 
produced photon polarization vectors, respectively, and a and b the eight-fold color indices 
of the off-shell gluons. Then the relevant matrix element can be presented as follows: 

9 / X ^ „ Pi — ki + nia \ k2 — P2 + h „ /V / , V 

- (pi - kiY - (fca - P2r 
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M2 = eg u{pi)f^ ey— -^7 e^— jr 7%'u(p2), (5 

9 / X ^ „ Pi — ki + ran ), „ —P2 — P + ran \ / x 

iTT'i - [Pi - hr i^q - {-P2 - pY 



M, = eg^u{p,)t'Ye. 'I /" ^ ^"7% /\ ' ^ "' '^ 'e,u{p^, (7) 

- (pi - k2Y '^m^g- {-p2 - pY 



9 / N X Pi + P + nia h „ ki — P2 + nio „ ,, , . 

"^g - (Pl + PY "^q - {kl - P2Y 

KA 2-r \ X Pi+P + ^q .a 1^ k - P2 + niq u ^ 

Ma = eg u{pi)-f^ex^ — 7 -^t Te^— 7 e^u{p2), (9) 



'm'^-{pi+pY ^rnl- {k2-p2, 



M7 = eg' u{p,) YC'^'iki, ^2, -ki - h) .'''l .. f'^fx 

{ki + k2Y 

-p2-p + mq ^ 

- y-P2 -pY 

KA 2-1 \ \ Pi+P + mq 

Ms = eg' u{pi) -f^ex— . ^ ' x 

fn'q - [Pi +pY 

xYc^'^'iki, k2, -k, - ^2) ,, Y;,, r''yt^(P2). 

(fci + k2Y 



(10) 



(11: 



In above expressions C'^^^{k,p, q) is related to the standard QCD three-gluon coupling 

C'^'iKp, q) = g^'iP - ky + /^(g - pY + g''{k qy. (12) 
The summation on the produced photon polarization is carried out by covariant formula 

J2e^e*'' ^ -g^r (13) 
In the case of initial off-shell gluon we use the BFKL prescription [11-13]: 

e^iky^ik)^^. (14) 

This formula converges to the usual one (13) after azimuthal angle averaging in the kx ^ 
limit. The evaluation of the traces in (4) — (11) was done using the algebraic manipulation 
system FORM [36]. The usual method of squaring of (4) — (11) results in enormously 
long output. This technical problem was solved by applying the method of orthogonal 
amplitudes [37]. 

The gauge invariance of matrix element is a subject of special attension in the kr- 
factorization approach. Strictly speaking, the diagrams shown in Fig. 2 are insufficient 
and have to be accompanied with the graphs involving direct gluon exchange between the 
protons (these protons are not shown in Fig. 2). These graphs are necessary to maintain the 
gauge invariance. However, they violate the factorization since they cannot be represented 
as a convolution of gluon-gluon fusion matrix element with unintegrated gluon density. The 
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solution pointed out in [10] refers to the fact that, within the particular gauge (13), the con- 
tribution from these unfactorizable diagrams vanish, and one has to only take into account 
the graphs depicted in Fig. 2. We have succesfuUy tested the gauge invariance of matrix 
element (4) — (11) numerically. 

2.3 Fragmentation contributions 

Perturbation theory becomes nonapplicable when the wavelength of the emitted photon 
(in the emitting quark rest frame) becomes larger that the typical hadronic scale O (1 
GeV~^). Then the nonperturbative effects of hadronization or fragmentation must be taken 
into account. Accordingly, the calculated cross section can be split into two pieses 

da = c?a-^ircct(/^^) + '^'^fragm(/^^) 

with dcr (XiYQcti^^'^) representing the perturbative contribution and da^^^^^[iji?) the fragmen- 
tation contribution. In our calculations we choose the fragmcntationn scale /i^ to be the 
invariant mass of the quark + photon subsystem, /x^ = {p + Pi)"^, and restrict c?cr(iirect(A*^) 
to /X > M ~ 1 GeV. Under this condition, the contribution (icrjjj.g^^(/i^) is free from diver- 
gences, so that the mass of the hght quark ruq can be safely sent to zero. The sensitivity of 
our results to the choice of M is reasonably soft, as we will discuss in Section 3. As far as 
the fragmentation contribution is concerned, its size is dramatically reduced by the photon 
isolation cuts (see below) . According to the estimates presented in Ref . [38] , the contribution 
from (icTfj-g^gj^ amounts to about 10% of the visible cross section. This value is smaller than 
the theoretical uncertainty in calculating the perturbative contribution d(^(]^\Yecii ^^'^ 
neglected in our analysis. 

2.4 Photon isolation cuts 

In order to reduce huge background from the secondary photons produced by the decays 
of tt" and t] mesons the isolation criterion is introduced in the experimental analyses. This 
criterion is the following. A photon is isolated if the amount of hadronic transverse energy 
deposited inside a cone with aperture R centered around the photon direction in the 
pseudo-rapidity and azimuthal angle plane, is smaller than some value E^^'^: 

(^had_^)2^^^had_^)2<^2_ U^J 

The both D0 and CDF collaborations take it! ~ 0.4 and E^^ ~ 1 GeV in the experi- 
ment [1-6]. Isolation not only reduces the background but also significantly reduces the 
fragmentation components connected with collinear photon radiation. It was shown that 
after applying the isolation cut (15) the contribution from the fragmentation subprocesses 
is strongly suppressed [38] . 
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2.5 Cross section for prompt photon hadroproduction 

In general case, according to fcr-factorization theorem the cross section of prompt photon 
hadroproduction can be written as a convolution 



271 , , 

A[X2,\<^2T-.^^ )dii2Tl^da{g g qq'j), 

X2 2,71 

where d-{g*g* — > qq'j) is the partonic cross section, ^(x, k|,,/i^) is the unintegrated gluon 
distribution in a proton and (pi and 02 are the azimuthal angles of the incoming gluons. The 
multiparticle phase space Ild^pi/2Ei5^^\j2p^^~J2p"^^) is parametrized in terms of transverse 
momenta, rapidities and azimuthal angles: 

% = ^dp^,%^. (17) 

Using the expressions (16) and (17) we can easily obtain the master formula: 



25Q7T^{xiX2sy 

..A/ 12 2 s A/ ,2 2n n 2 J1 2 . 2 2 r r . #1 #2 #1 #2 

Ztt ztt ztt ztt 



where \Ai{g*g* i?^7)P is the off-mass shell matrix element squared and averaged over 
initial gluon polarizations and colors, ipi and ip2 are the azimuthal angles of the final state 
quark and antiquark, respectively. We would like to point out again that \A4{g*g* i?^7)P 
is strongly depends on the non-zero transverse momenta kfj, and k'^j'. If we average the 
expression (18) over 0i and 02 and take the limit k^^ — >■ and k|y 0, then we recover the 
expression for the prompt photon hadroproduction in the usual coUinear approximation. 

The multidimensional integration in (18) has been performed by the means of Monte 
Carlo technique, using the routine Vegas [39]. The full C-I--I- code is available from the 
authors on request^. 



3 Numerical results 

3.1 Theoretical uncertainties 

Except the the unintegrated gluon distribution in a proton A{x, k'^, /i^), there are several 
parameters which determined the overall normalization factor of the cross section (18): the 
quark mass m^, the factorization and normalisation scales and fi^. 

Concerning the unintegrated gluon densities in a proton, we have tried here two different 
sets of them. These sets are widely discussed in the literature (see, for example, review [35] 
for more information). Here we only shortly discuss their characteristic properties. 

^lipatovQtheory.sinp. msu.ru 



7 



First unintegrated gluon density used has been obtained [40] recently from the numerical 
solution of the full CCFM equation. Function .4(a;, k|,, /x^) is determined by a convolution 
of the non-perturbative starting distribution Ao{x) and the CCFM evolution denoted by 



In the perturbative evolution the gluon splitting function Pgg{z) including non-singular terms 
(as described in detail in [41]) is applied. The input parameters in Aq{x) were fitted to de- 
scribe the proton structure function ^2(2;, Q^). A acceptable fit to the measured F2 values 
was obtained with /ndf — 1.83 using statistical and uncorrelated systematical uncertain- 
ties (compare to 'x^ /ndf ~ 1.5 in the collinear approach at NLO). This distribution has been 
applied recently in the analysis of the deep inelastic proton structure functions F|, F| and 



Another set (the so-called KMR distribution) is the one which was originally proposed 
in [31]. The KMR approach is the formalism to construct unintegrated gluon distribution 
from the known conventional parton (quark and gluon) densities. It accounts for the angular- 
ordering (which comes from the coherence effects in gluon emission) as well as the main 
part of the collinear higher-order QCD corrections. The key observation here is that the fi 
dependence of the unintegrated parton distribution enters at the last step of the evolution. In 
the numerical calculations we have used the standard GRV (LO) parametrizations [43] of the 
collinear quark and gluon densities. The KMR-constructcd parton distributions were used, 
in particular, to describe the prompt photon photoproduction at HERA [44] and prompt 
photon hadroproduction Tevatron [32]. 

Significant theoretical uncertainties are connected with the choice of the factorization and 
renormahzation scales. The first of them is related to the evolution of the gluon distributions, 
the other is responsible for the strong coupling constant as{n'ji}- often done, we choose 

the renormahzation and factorization scales to be equal: = /.ip = fJ' = ^\pt\- In order to 
investigate the scale dependence of our results we will vary the scale parameter ^ between 
1/2 and 2 about the default value ^ = 1. 

The masses of all light quarks were set to be equal to ruq — 4.5 MeV and the charmed 
quark mass was set to rric = 1.4 GeV. We have checked that uncertainties which come 
from these quantities are negligible in comparison to the uncertainties connected with the 
unintegrated gluon distributions. For completeness, we use LO formula for the strong cou- 
pling constant Q!s(/x^) with n/ = 3 active quark fiavours at Aqcd = 232 MeV, such that 
q;s(M|) = 0.117. Note that we use special choice Uf — A and Aqcd = 130 MeV in the case 
of CCFM gluon (^^(Ml) = 0.118), as it was originally proposed in [40]. 

3.2 Inclusive prompt photon production at Tevatron 

The experimental data [1-6] for the inclusive prompt photon hadroproduction — > 7 X 
at Tevatron come from both the D0 and CDF collaborations. The D0 [1, 2] data were 
obtained in the central and forward pseudo-rapidity regions for two different ccntcr-of-mass 
energies, namely ^/s = 630 GeV and ^/s = 1800 GeV. The central pseudo-rapidity region is 
defined by the requirement \r]\ < 0.9, and the forward one is defined by 1.6 < \r]\ < 2.5. The 
more recent CDF data [3, 4] refer to the same central kinematical region I77I < 0.9 for both 




(19) 



Fl [42]. 
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beam energies y/s = 630 GeV and y/s = 1800 GeV. As usually, the photon pseudo-rapidity rj 
is defined as 77 = — lntan(6'/2), where 6 is the polar angle of the prompt photon with respect 
to the proton beam. Also the CDF collaboration has presented a measurement [5] of the 
prompt photon cross section at ^/s = 1800 GeV which is based on events where the photon 
converts into an electron-positron pair in the material of inner detector, resulting in a two 
track event signature ("conversion" data). These data refer only to the central kinematical 
region. Actually, there are two different datasets, which were used in the CDF measurement 
with conversions, namely 8 GeV electron data and 23 GeV photon data^. Last available 
experimental data has been presented by the D0 collaboration very recently [6]. These data 
extend previous measurements [1-5] to significantly higher values of photon pt (namely to 
Pt ~ 300 GeV at = 1960 GeV). 

The results of our calculations are shown in Figs. 3 — 18. So, Figs. 3 and 4 confront 
the double differential cross sections da/dExdrj calculated at ^/s = 630 GeV in different 
kinematical regions with the D0 [2] and CDF [3] data. The solid and dottes lines are obtained 
by using the CCFM and KMR unintegrated gluon densities with the default factorization and 
renormalization scale. The upper and lower dashed lines correspond to the scale variations 
in CCFM gluon as it was described above. The data/theory ratio is also shown in Figs. 5 
and 6. One can see that our predictions agree well with the experimental data within the 
scale uncertainties. However, the results of calculation with the KMR gluon at the default 
scale tend to underestimate the data in the central kinematical region and agree with the 
D0 data in the forward rj region. The CCFM gluon density gives a perfect description of 
the data in both kinematical regions. Main difference in predictions of the KMR and CCFM 
gluon densities concentrated in the central pseudo-rapidity region at low pt- This region 
corresponds to the low x values where effects of the gluon evolution play the dominant 
role. It is clear from the Figs. 7 and 8 where gluon and valence quark contributions are 
shown separately. In the central pseudo-rapidity region the gluon-gluon fusion dominates up 
to Pt ~ 20 GeV approximately whereas in the forward pseudo-rapidity region the valence 
quark contribution should be taken into account everywhere. In these plots we have used the 
CCFM gluon density for illustration. Note that the KMR gluon give description of the data 
which is rather similar to collinear NLO QCD calculations [17]: the results of measurement 
are higher than the NLO prediction at low Et in the central rj range but agree at all Et in 
the forward pseudo-rapidity region. 

It is important that KMR predictions based on the g*g* —>■ qq'-f off-shell matrix element 
practically coincide with the previous ones [32] which were based on the quark-gluon fusion 
and quark- ant iquark annihilation. This fact gives an additional check of our calculations 
and the self-consistency of the whole KMR scheme. 

The double differential cross sections da/dETdrj compared with the experimental data 
at i/i = 1800 GeV in different pseudo-rapidity regions are shown in Figs. 9 — 11. The 
data/theory ratio is depicted in Figs. 12 and 13. The gluon and valence quark contributions 
to the cross section at this energy are shown separately in Fig. 14 and 15. We find that our 
predictions with CCFM gluon density reasonable agree with the D0 [1] and CDF [3-5] data 
both in normalization and shape. However, the using of the KMR-evolved gluon density 
results to underestimation of the measurement by a factor about 2 in the central r] region. 

^See Ref. [5] for more details. 
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The similar situation is observed at a/s = 1960 GeV (see Fig. 16). This is due to different 
small-x behaviour of the unintegrated gluon densities, since valence quark contribution in 
this region is about 7% only, as it was demonstrated in Fig. 14. In forward pseudo-rapidity 
region the results obtained with the both gluon densities under consideration agree well with 
the experimental data. So one can conclude that in the central pseudo-rapidity region at 
low pt the observed cross section is a strongly sensitive to the unintegrated gluon density 
used. Therefore available experimental data can be applied in future to better constraint 
of the unintegrated gluon distribution. Specially it will useful in the case of the CCFM 
evolution (which is valid for both small and large values of x) since measurements [1-6] on 
the prompt photon production refers to whole x range. One can see that the gluon-gluon 
contribution at \/s = 1800 GeV dominates up to pt ~ 50 GeV in the central t] region and 
up to Pt ~ 20 GeV in the forward one, where x is small still. In general, the contribution 
from valence quarks is more significant in the forward rapidity region. This situation is the 
similar to one at y/s = 630 GeV. 

Concerning the coUinear approximation of QCD, it was shown [6] that results from NLO 
calculations agree with the recent high px measurement within uncertainties. However, at 
moderate and small pt region (which corresponds to small values of x) the shape of the 
measured cross sections at y/s = 1800 GeV is generally steeper than that of the NLO 
predictions. It was claimed [3, 4] that this shape difference is difficult to explain simply by 
changing the renormalization/f actor ization scales in the collinear calculation, or the set of 
parton distribution functions. 

Now we can try to extrapolate our theoretical predictions to LHC energies. We perform 
the calculation for both central and forward photon pseudo-rapidities rj. As a representative 
example, we will define the central and forward kinematical regions by the requirements 
|?7| < 2.5 and 2.5 < |?7| < 4, respectively. The transverse energy Et distributions of inclusive 
prompt photon production at y/s = 14 TeV are shown in Figs. 17 and 18. The formulas (15) 
were used in these predictions. These Figs show the size of theoretical uncertainties connected 
with the unintegrated gluon densities. It is worth mentioning that the extrapolation of the 
available parton distribution to the region of lower x is a special problem at the LHC energies. 
In particular, one of the problem is connected with the correct treatment of saturation effects 
in small x regioij^. Therefore additional work needs to be done until these uncertainties will 
be reduced. 

4 Conclusions 

We have tried a new theoretical approach to the production of prompt photons in hadronic 
collisions at high energies. Our approach is based on the fcf-factorization scheme, which, un- 
hke many early calculations [21, 22, 25], provides solid theoretical grounds for adequately 
taking into account the effects of initial parton momentum. The central part of our consid- 
eration is the off-shell gluon-gluon fusion subprocess g*g*—>-qq'j- The corresponding off-shell 
matrix elements have been calculated for the first time. At the price of considering the 2— >-3 
rather than 2— >2 matrix elements, we have reduced the problem of unknown unintegrated 
quark distributions to the problem of gluon distributions. This way enables us with making 

^See, for example, Ref. [35] 
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comparisons between the different parton evolution schemes and parametrizations of parton 
densities, in contrast with previous calculations [18, 32] where such a comparison was not 
possible (for the lack of unintegrated quark distributions except KMR) . 

Since the gluons are only responsible for the appearance of sea, but not valence quarks, 
the contribution from the valence quarks has been calculated separately. Having in mind that 
the valence quarks are only important at large x, where the traditional DGLAP evolution is 
accurate and reliable, we have calculated this contribution within the usual coUinear scheme 
based on 2 — > 2 partonic subprocesses and on-shell parton densities. 

We have calculated the total and differential cross sections and have made comparisons 
to the recent D0 and CDF experimental data. In the numerical analysis we have used 
the unintegrated gluon densities which are obtained from the full CCFM equation as well as 
from the KMR prescription. We have found that in the central pseudo-rapidity region at low 
Pt the observed cross section is strongly sensitive to the unintegrated gluon density. In this 
kinematic region the cross sections evaluated with CCFM and KMR gluon distributions differ 
from each other by a factor of about 2. This is due to the fact that low-pr measurements 
in the central pseudo-rapidity region probe the small-x region, where the shapes of the 
considered gluon densities are different. In the forward pseudo-rapidity region the predictions 
from different gluon densities practically coincide. We have demonstrated that the available 
experimental data can be used to constrain the unintegrated gluon distributions. Especially, 
it may be useful in the case of CCFM evolution (which is valid for both small and large 
values of x) since the recent D0 and CDF measurements cover the whole x range. 

Our results based on Kimber-Martin-Ryskin [31] gluon density and off-shell gluon-gluon 
fusion matrix element g*g*^qq'j agree with the previous results [32] based on the quark- 
gluon fusion qg^q^ and quark- antiquark annihilation qq-^g^ matrix elements and KMR [31] 
off-shell quark and gluon distributions. This can be regarded as an additional proof of the 
consistency of the proposed method. 
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Figure 1: Kinematics of the g*g* — > qq^ process. 
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Figure 2: Feynman diagrams which describe the partonic subprocess g*g* — > qqj at the 
leading order in ag and a. 
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Figure 3: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at \r}\ < 0.9 and ^/s — 630 GeV. The solid hne corresponds to the CCFM 
gluon density with the default scale fi = Et, whereas upper and lower dashed lines corre- 
spond to the usual scale variation in the CCFM distribution, respectively. The dotted line 
corresponds to the KMR unintegrated gluon density. The experimental data are from D0 [2] 
and CDF [3]. 
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Figure 4: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at 1.6 < \r)\ < 2.5 and s/s = 630 GeV. Notations of all curves are the 
same as in Figure 3. The experimental data are from D0 [2]. 
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Figure 5: The data/theory ratio for inclusive prompt photon hadroproduction calculated at 
\r]\ < 0.9 and ^/s = 630 GeV. Black and open circles correspond to the CCFM and KMR 
gluon densities, respectively. The experimental data are from D0 [2] . 
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Figure 6: The data/theory ratio for inclusive prompt photon hadroproduction calculated at 
1.6 < \r]\ < 2.5 and s/s = 630 GeV. Notations are the same as in Figure 5. The experimental 
data are from D0 [2] . 
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Figure 7: The double differential cross section da / dExdri for inclusive prompt photon 
hadroproduction at I77I < 0.9 and \/s = 630 GeV. The dashed and dotted lines represent 
the gluon and valence quark contributions to the prompt photon cross section, respectively. 
The sohd line corresponds to the sum of these contributions. We have use here the CCFM 
gluon densities. The experimental data are from D0 [2] and CDF [3] . 
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Figure 8: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at 1.6 < \r)\ < 2.5 and s/s = 630 GeV. Notations of all curves are the 
same as in Figure 7. The experimental data are from D0 [2] . 
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Figure 9: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at \r)\ < 0.9 and = 1800 GeV. Notations of all curves are the same 
as in Figure 3. The experimental data are from D0 [1] and CDF [3, 4]. 
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Figure 10: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at < 0.9 and ^/s — 1800 GeV. Notations of all curves are the same as 
in Figure 3. The experimental data are from CDF [5]. 
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Figure 11: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at 1.6 < \r]\ < 2.5 and ^/s = 1800 GeV. Notations of all curves are the 
same as in Figure 3. The experimental data are from D0 [1]. 
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Figure 12: The data/theory ratio for inclusive prompt photon hadroproduction calculated 
at \r)\ < 0.9 and -^/s = 1800 GeV. Notations are the same as in Figure 5. The experimental 
data are from D0 [1] . 
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Figure 13: The data/theory ratio for inclusive prompt photon hadroproduction calculated 
at 1.6 < \ri\ < 2.5 and ^/s = 1800 GeV. Notations are the same as in Figure 5. The 
experimental data are from D0 [1]. 
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Figure 14: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at [t^I < 0.9 and ^/s = 1800 GeV. Notations of all curves are the same as 
in Figure 7. The experimental data are from D0 [1] and CDF [3, 4]. 
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Figure 15: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at 1.6 < \r]\ < 2.5 and ^/s = 1800 GeV. Notations of all curves are the 
same as in Figure 7. The experimental data are from D0 [1] . 
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Figure 16: The double differential cross section da/dETdr] for inclusive prompt photon 
hadroproduction at [t^I < 0.9 and -^/s = 1960 GeV. Notations of all curves are the same as 
in Figure 3. The experimental data are from D0 [6]. 
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Figure 17: The transverse enery Et distribution of inclusive prompt photon hadroproduction 
calculated at |7^| < 2.5 and ^/s — 14 TeV. Notations of all curves are the same as in Figure 3. 
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Figure 18: The transverse enery Et distribution of inclusive prompt photon hadroproduction 
calculated at 2.5 < < 4 and ^/s = 14 TeV. Notations of all curves are the same as in 
Figure 3. 
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